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Abstract 

The application of enhanced quasi-sparse eigenvector methods (EQSE) to the Hub- 
bard model is attempted. The ground state energy for the 4x4 Hubbard model is 
calculated with a relatively small set of basis vectors. The results agree to high 
precision with the exact answer. For the 8x8 case, exact answers are not available 



lO i but a simple first order correction to the quasi-sparse eigenvector (QSE) result is 

^ ■ presented. 

o 
o 
o 

^ ! 1 Introduction 

^ . 

^ ■ The enhanced quasi-sparse eigenvector (EQSE) method of solving quantum 

Q ! field theory Hamiltonians is the combination of quasi-sparse eigenvector (QSE) 

method [1] with a stochastic calculation for the contribution of the remaining 
basis states [2]. The Hubbard model was chosen as a laboratory for testing this 
;h ■ method for several reasons. First, we believed the approach could yield results. 

^ The basis vectors can be specified in a few words of data and the Hamiltonian 

is sparse is momentum space. On the other hand, the ground state of the 
Hubbard model is known for its inclusion of a extraordinary number of Fock 
states [3] so the model presents a non-trivial challenge to the quasi-sparse 
approach. Finally, the Hubbard model is thought to be a physically relevant 
model for superconductivity [3]. While the description of the model is simple, 
solutions have been difficult and any promising new approach is worthwhile. 

We work on a 2-dimensional spatial lattice with the Hubbard Hamiltonian 

< i,j > nearest 
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There are 2 species of electron so there are 4 possible states for each lattice site. 
Thus the 8x8 Hubbard model has 4^^^ dimensions. Even after using particle 
conservation to partition the space there arc more than 10^^ basis vectors. 
In this large space the Hamiltonian is clearly sparse But the equality of the 
off-diagonal elements contributes to an extraordinary number of Fock states 
in the ground. 

It is thought that the D-wave correlator Cd^2_y2 (^) important indicator 
of superconductivity [4] . 



2 Hamiltonian Momentum Lattice Formulation 



After making space periodic we use the Fock states as our basis set. The 
Hamiltonian conserves momentum which helps limit the number of relevant 
basis states. 

i^kin = -2t 2^(cos — + COS -^)<p,^a„p,^ 

k — l + m — n = 
p — q + r — s = 



There is a 16 member symmetry group generated by reflections in the x and 
y planes, x ^ y and|^|. For the purpose of finding the ground state, we use 
only symmetrized basis states. 

The first step in the calcualation consists of picking a basis set of size N and 
diagonalizing its submatrix of the Hamiltonian using the Lanczos method [5]. 
The ^basis vectors which least contribute to the ground are then discarded, 
replacements are chosen and the Hamiltonian is again diagonalized. When the 
ground energy Eq obtained in this manner converges the QSE step is complete. 

The EQSE is a first order correction to this result. We calculate it stochasti- 
cally. Let C be a set of basis vectors for the complement of the A^- dimensional 
subspace. Then, choosing representative basis vectors f e C with probability 
P{v) . 



Eq + Average — -V, FOiv) = -— — ^^-^ — — 

P[v) A„ 

where |0) is the ground state of the N dimensional subspace. The expectation 
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is thus 



FO{v) 



3 Results 



Results were obtained for the ground state energy, wavefunction and d-wave 
correlator. The computing time was about 2 days on a 350Mhz Pentium 11. 
Where available, we compare with Husslein et al [4] results labeled Exact, 
Projector Quantum Monte-Carlo (PQMC), and Stochastic Diagonalization 
(SD) (which uses a different method of choosing the subspace than QSE). 



Ground State Energy 4x4 Hubbard Model(^ occupied) 



Coupling 


States 


QSE 


EQSE 


Exact 


SD 


PQMC 




50 


-.47471 


-.50127(5) 








U=2 


100 


-.47967 


-.50147(5) 


-.50194 


-.501 


-.49 




500 


-.49454 


-.50181(3) 










1000 


-.50062 


-.50198 










50 


-1.5707 


-1.80635(5) 








U=4 


100 


-1.6203 


-1.8113(4) 


-1.8309 


-1.829 


-1.8(2) 




500 


-1.7476 


-1.8242(3) 










1000 


-1.8003 


-1.8302 










50 


-2.2450 


-2.6663(8) 








U=5 


100 


-2.3322 


-2.6724(7) 


-2.7245 


-2.723 


-2.9(3) 




500 


-2.5578 


-2.7073(4) 










1000 


-2.6512 


-2.7208(2) 










2000 


-2.685 


-2.7231 










50 


-2.963 


-3.615 










100 


-3.103 


-3.635 










1000 


-3.452 


-3.697 










2000 


-3.595 


-3.723 
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Ground Energy 8x8 Hubbard Model(|| occupied) 



Coupling 


States 


QSE 


EQSE 




50 


-.281 


-2.58(2) 




100 


-.443 


-2.49(2) 


U=2 


500 


-1.221 


-2.40(2) 




1000 


-1.751 


-2.406(8) 




2000 


-1.956 


-2.423(9) 




4000 


-1.958 


-2.427 




50 


-.811 


-9.18(6) 




100 


-1.386 


-8.48(6) 


U=4 


500 


-3.449 


-7.58(5) 




1000 


-4.798 


-7.59(3) 




2000 


-5.374 


-7.620(4) 




4000 


-5.387 


-7.621(5) 




50 


-1.222 


-13.33(8) 




100 


-1.88 


-12.33(8) 




500 


-4.65 


-10.65(6) 




1000 


-6.44 


-10.65(3) 




2000 


-7.225 


-10.671(5) 




4000 


-7.236 


-10.662(8) 




50 


-1.6 


-18.2 




100 


-2.4 


-16.6 




500 


-5.9 


-13.93(8) 




1000 


-8.15 


-13.78(4) 




2000 


-9.110 


-13.888(7) 




4000 


-9.113 


-13.880(8) 



The dx2_y2 correlator was also obtained using the QSE algorithm. The 4x4 
result again matched that of Husslein et al [4] . The EQSE calculation has not 
been completed and we therefore omit the data. 
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4 Conclusion 



As we can see, the ground state for the 4x4 Hubbard model can be well 
described with about 1000 symmetrized states and the 50 state results yield 
remarkable accuracy when the first order correction is included. In the 8x8 case 
it is clear that even 4000 states are not sufficient to describe the ground state. 
The precision of the first-order values will be determined with the completion 
of higher order calculations. 

Further advances will come in the refinement of the enhancement technique. 
Better importance sampling will speed convergence of second and higher orders 
contributions. Other extensions will be calculation of excited states of the 
Hamiltonian, correlation functions, binding energies and other quantities of 
interest. 

Acknowledgments The author thanks the organizers of the Light-Cone 
Meeting at Heidelberg, Nikolay Prokofyev and collaborators on the original 
papers cited here. 



References 

[1] D. J. Lee, N. Salwen, D. D. Lee, |hep-th/000225T . 

[2] D. Lee (in these proceedings) 

[3] P. W. Anderson, Science 235, 1196 (1987) 

[4] T. Husslein, W. Fettes, I. Morgenstern, |cond-mat /9705026| 

[5] Using Arpack package (www.caam.rice.edu/software/ARPACK). 



5 



